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We discuss the behavior of free perfect staggered fermions and truncated versions 
thereof. The study includes flavor non-degenerate masses. We suggest a new blocking 
scheme, which provides excellent locality of the perfect lattice action. A truncation proce- 
dure adequate for the structure of staggered fermions is applied. We consider spectral and 
thermodynamic properties and compare truncated perfect actions, Symanzik improved 
and standard staggered fermions in two and four dimensions. 
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1 Introduction 



The inclusion of dynamical fermions in simulations of gauge theories, in particular QCD, 
is a major issue in lattice field theory. Such simulations have generally been performed 
either using Wilson fermions |l| or staggered fermions ||. The latter formulation is 
especially useful in the chiral limit, because the remnant chiral symmetry U(l) ® U(l) 
protects the zero fermion mass from renormalization. As a related virtue, its artifacts due 
to the lattice spacing a are only of 0(a 2 ), whereas they are of 0(a) for Wilson fermions 
interacting by gauge fields. 

It is now widely accepted that the above lattice actions should be improved, so that 
the lattice spacing artifacts are suppressed and coarser lattices can be used, particularly in 
QCD simulations f3[. There are essentially two improvement strategies in the literature. 
In Symanzik's program M the action is improved order by order in a. For QCD with 
Wilson fermions this has been realized on-shell to the first order on the classical level || 
and recently also on the quantum level [Q. Less work has been devoted to the improve- 
ment of staggered fermions, perhaps because the artifacts in the standard formulation are 
already smaller. However, S. Naik has applied Symanzik's program on-shell, where he 
improved the free staggered fermion by adding more couplings along the axes 0, and the 
Bielefeld group did the same by adding diagonal couplings ||. Furthermore, the MILC 
collaboration achieved a reduced pion mass by treating the gauge variable as a "fat link" 



Finally, some work on improved operators has been done in this framework [10 



The other promising improvement scheme is non-perturbative in a and uses renormal- 
ization group concepts. It has been known for a long time that there are perfect lattice 
actions in the parameter space, i.e. actions without any lattice artifacts JITJ. Recently 



it has been suggested to approximate them for asymptotically free theories as "classi- 
cally perfect actions" [|12[| , which works very well in a number of two dimensional models 

For free or perturbatively interacting fields, perfect actions can be constructed ana- 
lytically in momentum space. For Wilson type fermions this has been carried out to the 
first order in the gauge coupling in the Schwinger model |l5j and in QCD [|l(J. A tech- 
nique called "blocking from the continuum" was extremely useful for this purpose. One 
expresses all quantities in lattice units after the blocking, and sends the blocking factor 
to infinity. Hence the blocking process starts from a continuum theory, and it does not 
need to be iterated in order to identify a perfect action. 

For staggered fermions a block variable renormalization group transformation (RGT), 
which does not mix the pseudo flavors and which does therefore preserve the important 
symmetries, has been suggested in Ref. |17]]. It requires an odd blocking factor n. Iterating 
the n = 3 block variable RGT, a fixed point action, i.e. a perfect action at infinite 



correlation length, has been constructed |13|, [L8|. Also for staggered fermions blocking 
from the continuum is applicable jnj. This has been carried out for a general (flavor 
non-degenerate) mass term, revealing the intimate relation to the Dirac-Kahler fermion 
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formulation in the continuum [p0|1 , and also including a suitable treatment of the gauge 
field ||21[| . Using the generalization to a flavor non-degenerate mass, the spectral doublers 
inherent to the staggered fermion formulation might be treated as physical flavors in a 
QCD simulation. 

For any blocking scheme, the perfect action displays the continuum scaling behavior by 
definition. In practice, however, the approximations in the determination of an applicable 
quasi perfect action violate this property to some extent. With this respect, the choice of 
the RGT is very important. 

In the present paper, we first review the procedure of blocking staggered fermions from 
the continuum, adding new aspects. We then discuss the optimization of locality, in the 
sense of an extremely fast exponential decay of the couplings in coordinate space. This 
property is crucial for practical applications, because we have to truncate the couplings to 
a small number, which is manageable in simulations. For excellent locality, the violation 
of perfectness due to truncation is not too harmful. Moreover, the virtues of staggered 
fermions come into play: the truncation induces errors of 0(a 2 ) at most, and the remnant 
axial symmetry is preserved in the truncated perfect action, hence it still excludes additive 
mass renormalization. On the other hand, a strong mass renormalization is a severe 
problem for truncated perfect fermions of the Wilson type [22], 123 . 



We suggest a new blocking scheme, which leads to a higher degree of locality than the 
usual block average method. Still the couplings decay exponentially in coordinate space. 
We then truncate them to a short range by means of mixed periodic and antiperiodic 
boundary conditions, which are particularly suitable for staggered fermions. We finally 
compare the spectra and the scaling of thermodynamic quantities for several truncated 
perfect, Symanzik improved and standard staggered fermions. The discussion of the 
spectrum includes also the case of non-degenerate masses of the pseudo flavors. It turns 
out that the new blocking scheme, which we call "partial decimation", does consistently 
yield excellent results. Moreover, for that scheme, optimal locality is reached by the 
simple 5 function RGT in the massless case, in contrast to the block average scheme. 

Simulation results for the Schwinger model based on this action will be presented in 
a subsequent paper [pi . 



2 The Staggered Blocking Scheme 

The construction of perfect actions for free staggered fermions by blocking from the con- 
tinuum has been described in Refs. [ 2D , [HJ]. Let us briefly review this procedure. The 
staggered blockspin fermions are defined in two steps. First we transform the Nf = 2 d ^ 2 
flavors of continuum Dirac spinors ip b a {x) (a: spinor index, b: flavor index) into the Dirac- 
Kahler (DK) representation by f(x,H)^ These functions are considered as component 

■"■For the relation of the DK formulation of continuum fermions J2q| with staggered lattice fermions we 
refer to Ref . [Em . A relation to the present block spin transformation is discussed in Ref . |2(J . 
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functions of inhomogeneous differential forms 

$ = ^ v (x,H)dx H > dx " = dx^ 1 A . . . A dx» h , (2.1) 



H 



where H = {/xi, . . . ,fih}, //i < . . . < \ih is a multi-index. Transformation and inverse 
transformation read 

= 4=1^**)' ^ = 7 W 7" 2 --- 7 Mfe , (2-2) 
V40«0 = 4^72 (2-3) 

Second, we introduce a coarse lattice of unit spacing f = {y \ y^ = n^}, which is a 
sublattice of T = {y \ = n^/2}, with n M , G 5Z. The fine lattice points y are uniquely 
decomposed as (/2 is the unit vector in /i-direction) 

y = V + e H /2 , e H = E £ • (2-4) 

Thus the multi-index H(y) defines the position of a fine lattice point y with respect to 
the coarse lattice T. Now the blockspin variables <p(y) can be defined as averages of the 
component functions ip(x, H(y)), with a normalized weight U(x — y), JdxU(x — y) = 1, 
which is assumed to be even and peaked around x = y, 

4>{y) = -j=H ll {y) * J dx U(x - y) ^(x) . (2.5) 



"f ab 

This scheme has been proposed first in Ref. |19[| . Its peculiarity is that the staggered block 
centers depend on the multi-index H of the Dirac-Kahler component functions ip(x,H). 
Block average (BA) means in this case average over the overlapping lattice hypercubes 
[y] = {x\ — 1/2 < (x M — y^) < 1/2}. This scheme is given by II = Hba, ^-ba(x) = 1 for 
x G [y], Uba(x) = otherwise. 

For the following calculation we diagonalize the lattice action using the staggered sym- 
metries. Here it is important that fine lattice shifts are no symmetry transformations. 
However, combination with site-dependent sign factors gives rise to the non-commuting 
flavor symmetry transformations |27]]. Therefore we replace ordinary Fourier transfor- 
mation by harmonic analysis with respect to flavor transformations and coarse lattice 
translations. We thus obtain a modified momentum representation which intertwines 
Fourier transformation and the transition back from DK fermions to the Dirac basis, 

€ip) = E ^7^0(2/), 

y 

«(p) = E ^7^0(2/), P eB=]-n,n] d . (2.6) 

y 
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Inserting Eq. (|2.5|) we find 



m = -i/^E W) €(p') E z l(p - pl)y ia b {y) T?? r , (2.7) 



where i^ b (p), n(p) denote the Fourier transform of ^(x), Il(x). The last sum can be 
re-written as 

E ^ E^ /2 t£ 7 a V = (2n) d E % - 2tt0 E II (- 1 )'" ^ 7 a V 
yer K ie% d K 

= N f (2n) d E S(q - 2*1) 7 fj° lT ^ ■ (2.8) 

We have used the orthogonality of the 7-matrix elements 

E7 a 1 1%* = N f 5 aa ,5 w , (2.9) 

H 

and H(l) is defined by H(l) = {/i|Z M is odd}, H = H for h even, if = {/i | p if} 
for fa odd. Finally Eq. (|2~7| ) becomes (summation over double spin and flavor indices is 
understood) 

ttto) = y/N~f E n(p+27r0 ^:(p+2tt0 , i> b (p+2nl) = 7 fj° ^{p^l) 7 *« t . (2.10) 
ie% d 

Note that the blockspin transformation is diagonal with respect to spin and flavor for 
continuum momenta within the first Brillouin zone B, yet not for all / 7^ 0. 

We are now prepared to compute the perfect action for a RGT of the Gaussian type. 
Starting from a continuum action with a general mass term m& - which does not need to 
be flavor degenerate - the perfect lattice action S[(f>, 0] is defined as 

e -sm = Jvjtyjvfp*, ex P {-|^-^^(- g )(i7^^ + m 6 )^(?)} 

X 6XP {/ ~N%V I m ~ P) " E n(p+27r/) $ a b (-p-2nl) ] V b (p) 



+V b a (-P) [4> b a (p) - y/Nf E n(p+2vr/) ^ a b (p+2nl) ] 

l£% d 

+V b a(-p)D b a Ap)V b a>(p)}} (2-11) 

A non-zero term 

DUp) = lL'Dl{p) + 5 aa >D b (p) (2.12) 
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"smears out" the blockspin transformation, as in Ref . [[Hj] . This term is used to optimize 
locality of the resulting perfect action; it will be specified later on. The Gaussian integrals 
over ■0,-0 and 77,77 can be evaluated by substitution of the classical fields leading to f\ 

S[M\ = I f tt(-P) G~Z'(P) €(P) 1 , (2-13) 



b N f (2tt) 



with the lattice propagator 



t (2 ' 14) 

Note that G is flavor diagonal, G 6b = G b <5 feb , because ^y H ")^ K j H ( l > is diagonal iff 7^ is. 
In particular, for a degenerate mass term the adjungation with j H w is trivial, and the 
lattice propagator is proportional to 8w in flavor space. We define 



GUP) = -*Y.laa>Q» + 5 ««'Q b a(P) > ( 2 ' 15 ) 

hence Q b Q become 

Q» = Blip) + E n(p+2 7 r0 2 ^EE^(0Tg7^/ ( f I^ a + ^ ' (2 ' 16) 

Q b o(p) = D b ( P ) + E n(p+27r0 2 ^EE^(07^^/ r . , , ■ (2.17) 

In case of a non-degenerate mass the sums J2v J2k can n °t be contracted according 
to Eq. (|2.9|) , due to the sign factor sk{1) = Y\u^k{~ !)'"■ However, in the degenerate case 
nif, = m, with flavor independent smearing terms D b = D, we simply obtain 

QM - dm + £ n^)' ffffiffi . (2.1s) 

«.&>> = + ("SI 

The perfect action in real space arises from Eq. Q2.13| ) inserting the momentum rep- 
resentation Eq. (|2.6|) . After some 7-matrix algebra 1P0|| , we arrive at 

S$,<f>) = Y,Hy)™{y,y')<f>(y') , (2.20) 
m(y,»0 = 'E I P K (1/)P(V-1/>1/) MK (V-1/)- ( 2 - 21 ) 

A" 



2 We ignore constant factors in the partition function. 
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Corresponding to a lattice propagator diagonal in flavor space, the sum over K runs over 
multi-indices K e V with diagonal j K , in the Weyl basis V = {0, 1 2, 3 4, 1 2 3 4}. The 
sign factors p(z,y) = p(H(z),H(y)) arise from ^ H ^ K = p(H, K) ■y HAK , where HAK = 
(H U K)\(H n if), and p K {y) is given by p(H(y), K)p(K, H {%/)). By symmetry, the only 
non-zero contributions to M K (y) are 

M* (y) = ipfo if) jf J^j e~*y (p) for ff(y) = pAJT , (2.22) 
M K (y) = jf * e HW M ir(p) for H(y) = K, (2.23) 

^ M» = (, 24) 

The flavor degenerate case leads to vanishing components for K ^ 0, and for M® = M^o 
we simply obtain 



It has been proven in Ref. []20| that the couplings given by the fermion matrix m(y, y') 



are local, i.e. they decay faster than any power of \y — y'\. For that, certain periodicity 
properties apply, which translate into (for simplicity of notion let p denote either p or 0, 
if AO = K) 

QiM = E •£ OfM . Offr + = { ^ [^^^ . (2.26) 

Again, we sum over diagonal 7-matrices only. It is provided that the corresponding re- 
quirements are met for the smearing terms D b ~^ within Qj~(p), see below. In consequence, 
the fermion matrix components M^(p) obey periodicity conditions analogous to Qa(p) 



/i \r / ^^^J I'^^w^^v.^vj ~^*.u.*^v^v^u.^ v^u-v^v^v^v*.^ v,v^ "X, fj, 

and the integrands of Eqs. ( |2.22| , p . 2 3|) are periodic with respect to the Brillouin zone B 



and analytic in a strip around the real axis. This implies locality of the perfect action. 

The coupling of even and odd lattice points is due to the components of the 
fermion matrix; the Mq components couple even-even and odd-odd. We add without 
proof that even-odd decoupling of the Hermitian matrix m can be shown in any even 
dimension d with arbitrary (non-degenerate) mass terms for truncated versions of the 
perfect fermion matrix m(y,y'), see Ref. |2(J for d = 2. This is a useful property in 
simulations with Hybrid Monte Carlo algorithms. However, with (non-perfect) coupling 
to a gauge field and non-zero mass term (i.e. with even-even, odd-odd as well as even-odd 
couplings), this is not true in general. Yet, since even-odd decoupling of rrv m is a perfect 
property, it might be imposed as a construction requirement for a (quasi) perfectly gauged 
fermion vertex. 
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3 Optimization of Locality 



In the blocking scheme described so far, there is quite some freedom left. In particular, 
we may use averaging functions different from IT^, and we can choose the smearing term 
D in Eq. ( |2.11| ). In both cases we aim at optimization of the locality in the resulting 



perfect action. 

Let us first discuss the averaging scheme. In Eq. ( |2.11| ) we implicitly assumed the 
same blocking of ip and ip, given by the weight function IT(a;) resp. its Fourier transform 
Il(p). Now we consider the case of a block average II = Hba for ip (?/>) only, while ip (i/j) is 
put on the lattice by decimation, IL(x) = 5(x). Thus we obtain a single factor IT(p+27rZ) 
in Eq. (fZTJ) with 

n(p) = u BA ( P ) = n - . h = 2sin T • (3- 1 ) 

In case of a 5-blockspin transformation (D = 0), this means to identify the averaged 
continuum and lattice 2-point functions 

( <Kv) 4>W) ) = I dx{<p(x,H(y))p(y',H(y'))) . (3.2) 
J [y] 

Due to translation invariance it doesn't matter whether we average source or sink of the 
continuum expression, or whether we allocate the space directions to be integrated over 
to source and sink in some way. The last point of view may be used to make a closer 
contact to the construction of staggered fermions from DK fermions in the continuum 



[26[, as discussed in Ref. p0[ . We call this blocking scheme partial decimation. For the 
2-point functions every space direction is integrated over once; therefore we do not run 
into the difficulties arising for blockspin transformations with complete decimation, which 
do not have a corresponding perfect action. 

For both blocking schemes, block average for ip and ip (BA) and partial decimation 
(PD), we now want to optimize locality of the couplings by making use of the smearing 
terms in Eqs. ( |2.16| , j2.17|) . As an optimization criterion it has been suggested to require 
that in the effectively Id case - with momenta p = (pi, 0, . . . , 0) - the couplings are 
restricted to nearest neighbors as in the standard action |n| 0-0 In the degenerate case 
we require 

i T iMi(pi,0,...,0) + Af (pi, 0, . . . , 0) = /(m)[ipi7i + to], (3.3) 
with m| m= o = and /(0) = 1. Our ansatz for the Gaussian smearing term reads 

D^p) = c{m)p^ D (p) = a{m). (3.4) 



3 It has been shown in Ref. j28] that this criterion does optimize locality in d = 4 for scalar particles 
over a wide range of masses. 
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Requirement (|3.3|) can be fulfilled in both blocking schemes we are considering, if we 
specify the RGT as follows 

Cba(^) = [cosh(m/2) — 1)] / m 2 , clba{™) = [sinh(m) — m] / m 2 , 
cp£)(m) = , ap£)(m) = [cosh(m/2) — 1] / m . (3.5) 

In both cases, we obtain rh = 2sinh(m/2). For m = a non-vanishing static smearing 
term a(0) would explicitly break the remnant chiral symmetry in the fixed point action. 
Therefore, the static term should vanish for optimized locality, as it does in both cases. 
Furthermore, in the PD scheme the chiral limit is optimized for locality by a simple 5 
function RGT. As an advantage of this property - which is not provided by the BA scheme 
- there is a direct relation between the n-point functions in the continuum and on the 
lattice. In addition, the extension to interacting theories might involve numerical RGT 
steps in the classical limit, which also simplify in the absence of a Gaussian smearing 
term. Finally, a non-vanishing term c(m) causes complications in the inclusion of a gauge 
interaction; the PD scheme avoids such problems. 

The decay of the couplings m(x, 0) = m(x) in the massless case for d = 2 and d = 4 
is shown in Figures 1 and 2. We see that the PD blocking scheme works better. In 



n - 



X 
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Figure 1: The decay of couplings in (1,0) direction for d = 2 in the massless case: BA 
without smearing (diamonds), BA optimized (squares), PD (triangles). 



the non-degenerate case the simplest ansatz for a smearing term is (corresponding to the 



staggered fermion action with non-degenerate mass [[29 

£>» = cp, , D b (p) = £ 7 * a K J] cos(p M /2) . 

Kev fiei< 



(3.6) 
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Figure 2: The decay of 4d couplings in various directions in the massless case: BA with 
optimal smearing term and PD. We see that the latter decays faster. 



In this case, we obtained - by means of numerical optimization - a similar decay of 
couplings as with degenerate masses, yet no strict 1-dimensional ultralocality. For d = 2 
it appears that a non-degenerate parameter a±2 ^ does not improve the coupling decay 
significantly, so we worked with an = 0. Again the PD scheme, where we assumed 
cpd = 0, led to a more local action, see Figure 3. 



4 Truncation Effects 

The litmus test of any blocking scheme is given by the truncation effects for a practicable 
number of remaining couplings. These effects are minimized by maximal locality. Yet, 
the truncation scheme itself may have some impact on the truncation errors. An elegant 
procedure has been proposed in Ref. for Wilson fermions. First the perfect action 
is constructed on a small lattice volume N d by restricting the momentum components 
to the discrete values = 2irn fl /N e] — 7r, it), 6 Z. Typically, one chooses N = 3, 
and the resulting couplings are then used also in a large volume too, where they are not 
exactly perfect any more. This truncation scheme has the virtues of automatically correct 
normalizations, and a simplification of the numerical evaluation as opposed to a truncation 
in coordinate space. Furthermore, the mapping to the corresponding truncated perfect 
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Figure 3: The decay of couplings in (1,0) direction for d = 2 with masses mi = 0, m2 = 1: 
BA without smearing (diamonds), BA optimized (squares), PD (triangles). 



action in a lower dimension remains exact; for instance, one can reproduce the effectively 
Id nearest neighbor action starting from d > 1 by summing over the extra dimensions, 
which provides a sensitive test of the numerical accuracy. Provided a good locality, this 
method works well for Wilson-type fermions, gauge fields JZ2J and scalars |[28|| . 

For staggered fermions, only a crude truncation in coordinate space has been applied 
so far [21]. One includes the coupling distances ±1/2, ±3/2 in the \i direction, and 



1,0,-1 in the non-/z directions. This involves more couplings than the N = 3 truncated 
Wilson-type fermion (called "hypercube fermion"), whereas the number of degrees of 
freedom is the same in both cases. However, although significant improvement has been 
achieved compared to the standard staggered formulation, the quality of its spectral and 
thermodynamic properties did not reach the level of the "hypercube fermion". In order 
to arrive at results of similar quality, we adapt the above truncation procedure - together 
with the PD scheme - to the case of staggered fermions. 

We treat the components of the fermion matrix M~f(p),jl = fi, separately, as they 
show different behavior under reflections, 

M K (n n n n n \ - / ~ M f(P) for V G K Afx 

M A ( Pl . . . Pv _ u - Pu ,p v+l ...p d ) - | M ^ (p) for v ^ RA ~ ■ (4-1) 

Again we use the short-hand notation KAjl = K for jl = 0. Truncation is achieved by 
discrete Fourier transformation and a discrete support given by cn(v^) = 1,1/2,0 for 
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|2t/ m | <, =, > N, respectively, 

MUV) = n^r^ 1 E e- iPy Mf{ P ). (4.2) 

It is the set of discrete momenta B^. N which depends on the reflection properties given 
by /2, K. We choose partially antiperiodic boundary conditions in y-space, 

D F b k ^ _ _ f (2n + 1)tt/JV for z/ G if A/2 „ , , , 

p G o^.jv ^ p„ - j 2n7r/jV for v ^ KAj2 , n G ^, p u Gj - 7r,7rj . (4.3) 



It is easily verified that the transformation of the truncated components back to momen- 

>p,;N > 



turn space reproduces the perfect values at p G Bj^.j 



M? N {p)=Y,j py M? N {y) = Mf{p) for p e B« N . (4.4) 

The components M~f(y) inherit the reflection behavior of Eq. ( |4.1| ) by Fourier transfor- 
mation. Therefore, treating them as periodic functions in y-space by discrete momenta 
p v = 2nn/N, n G 2l> would make them vanish at the boundaries y v = ±N/2 for v G KAfi 
artificially. We avoid this effect by the above choice of discrete momenta. It pays off by 
a drastic reduction of truncation effects, see Figure § for the 2d spectrum in the massless 
case using partial decimation. The values of the perfect couplings, truncated by mixed 
periodic boundary conditions as described above, are given in Table [I]. 

Furthermore, the truncation effects are significantly stronger for even truncation dis- 
tances N; in particular, for N = 2 and N = A the energies become complex-valued at 
large momenta. 

For iV = 3 - which we consider to be tractable in numerical simulations - the amount 
of improvement is already striking. We compare the spectra for d = 2 in the PD and BA 
blocking scheme in the massless case (Figure |j) and for non-degenerate masses mi = 
0,m 2 = 1 (Figure ^]). The couplings derived from the PD blocking scheme turn out to 
be better, in agreement with the higher degree of locality observed in the (untruncated) 
perfect action. We see that the improvement is still good in the case of non-degenerate 
masses. In this case, we optimized the smearing parameters numerically, as discussed 
in Section ||. The values are a B A = 0.060, cba = 0.125, apn = 0.041, cpo = 0. The 
standard staggered fermion results for non-degenerate masses are calculated with the 
action proposed in Ref. [p9|1 . We emphasize that the spectra of (untruncated) perfect 
actions are indeed perfect, i.e. identical to the continuum spectra (up to the periodicity, 
which is inevitable on the lattice). Hence the spectrum reveals directly the artifacts due 
to truncation. 

Figure 6 compares the massless spectra in d = 4, and we see that the qualitative 
behavior observed in d = 2 persists. We compare the standard staggered fermion, and the 
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Pi 

Figure 4: The massless d = 2 spectrum. We compare standard staggered (fat dots) and 
perfect spectrum (full line) with perfect actions truncated with N = 3 and mixed periodic 
boundary conditions: block average (dotted), partial decimation (dashed). The last case 
is also plotted for truncation with periodic boundary conditions (dashed-dotted) . 



optimized BA and PD fixed point fermions, both truncated by mixed periodic boundary 
conditions. We see again that the PD scheme is superior. Furthermore, we show for 
comparison also the dispersion relation of a Symanzik improved action called "p6" from 
Ref. M. Symanzik improvement by additional couplings along the axes (Naik fermion, [|7|]) 



only yields a moderate quality |^TJ , but the Bielefeld group suggests a number of actions, 
where Symanzik improvement is achieved by diagonal couplings. The p6 action is the 
best variant among them, and also Figure 6 confirms its excellent level of improvement. 
However, it is not obvious how this formulation can include a general mass term in a 
subtle way. 

As a further test in d = 4, we consider the thermodynamic ratio P/T 4 (P: pressure, 
T: temperature). According to the Stefan-Boltzmann law, this ratio is 77r 2 /180 for mass- 
less fermions in the continuum, and a lattice action with many discrete points N t in the 
temporal direction will asymptotically reproduce that value. However, the speed of con- 
vergence, and in particular the behavior at small Nt, depend on the quality of the action. 
In contrast to the spectrum, this ratio is not even exact for the fixed point action in 
Eqs. ([2.18 , 2.19 ), because of the "constant factors" that we ignored when performing the 
functional integral in the RGT Eq. fl2.11|) . Such factors may depend on the temperature, 



so with respect to thermodynamics our action is not fully renormalized ||22|| . However, it 
turns out that the unknown factor is very close to 1, except for the regime of very small 
Nt (about Nt < 3), which corresponds to extremely high temperature. So the main issue 
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Figure 5: The d = 2 spectrum for m\ = 0,m 2 = 1. Standard staggered (fat dots) 
and perfect spectrum (full line) compared with perfect actions truncated with N = 3 and 
mixed periodic boundary conditions: block average (dotted), partial decimation (dashed). 



is again the contamination due to the truncation. 

In Figure 7 we compare this thermodynamic scaling for a variety of staggered fermion 
actions at m — 0, and again the PD scheme turns out to be very successful. A similar 
level of improvement can be observed for the p6 action, which is also here by far better 
than the Naik fermion. We conclude that a good improvement method should in any 
case include diagonal couplings, which are far more promising than additional couplings 
on the axes (just consider rotational invariance, for example). Thus we can expand the 
range, in which a practically accurate continuum behavior is observed, by about a factor 
of 4 compared to standard staggered fermions (see Figures 4, 5, 6 and 7). 



5 Conclusions 

Spectral and thermodynamic results for truncated perfect staggered fermions have been 
presented before pl |, but the improvement did not really reach a satisfactory level there. 



We could now push that improvement significantly further, mainly thanks to the new 
blocking scheme, which we call partial decimation, but also with the help of a new trunca- 
tion technique {mixed periodic boundary conditions). We now reached a level of excellent 
improvement, similar to the results for truncated perfect Wilson-type fermions. 

As a novelty, we extended the construction of perfect actions to non-degenerate flavors, 
and we could preserve the same level of improvement after truncation also in that case. 
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a — 4 






a — z 








rU 


DA 




rU 


dA 


(0.5,0,0,0) 


0.696388 


0.663116 


(0.5,0) 


0.878234 


0.866299 


(0.5,1,0,0) 


0.040981 


0.045925 


(0.5,1) 


0.060883 


0.066850 


( U. 0,1,1, U ) 


n nn/i /i si 

U.UU4481 


n nn/i&7i 

U.UU4o 1 1 








(0.5,1,1,1) 


0.000495 


0.000361 








(1.5,0,0,0) 


0.015219 


0.023441 


(1.5,0) 


0.010425 


0.016043 


(1.5,1,0,0) 


-0.000431 


-0.000639 


(0.5,1) 


-0.005212 


-0.008022 


(1.5,1,1,0) 


-0.000767 


-0.001211 








(1.5,1,1,1) 


-0.000428 


-0.000635 









Table 1: The couplings of the perfect action for massless staggered fermions, constructed 
from the "partial decimation" (most successful) and from the optimized "block average" 
scheme, and truncated by mixed periodic boundary conditions. The couplings are odd in 
the half- integer component, and even in all other components. Among the latter there is 
also permutation symmetry. 



This is potentially important for the study of the decoupling of heavy flavors, or for QCD 
simulations with realistic quark masses. 

At m = 0, which is well described by staggered fermions, the PD scheme is optimally 
local if we just use a 6 function RGT. This simplifies the relation to the continuum n-point 
functions, and numerical RGT steps, which could be performed for interacting theories 
(in the classical limit). 

The next step is the inclusion of gauge interactions, which is not straightforward. 
However, also with this respect the absence of a Gaussian smearing term in the chiral 
limit is of advantage, because a momentum dependent smearing term - as present for the 
block average scheme - causes complications in that step. 

As a simple ansatz, one can insert standard link variables, connect the coupled sites 
over the shortest lattice paths and take the mean value over these paths. This method is 
a simplification, which is not perfect, but for Wilson type fermions it leads for instance to 



a drastically improved mesonic dispersion relation ||22j| . In that framework, applications 



to the charmonium spectrum [ 30(1 , pionic systems [23], and the suitable generalization of 



preconditioning techniques for the fermion matrix |31|] are under investigation. 

We have implemented this method for truncated perfect staggered fermions in d = 2. 
Including also "fat links" , we obtained promising results for the scaling properties in the 
Schwinger model |M| down to small inverse couplings, /3 < 1. In particular, the "pion 
mass" is drastically reduced compared to the staggered standard action, and the "eta 
mass" follows the prediction by asymptotic scaling very closely. 

Acknowledgment We thank S. Chandrasekharan for his assistance. 
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continuum 




Figure 6: The dispersion relation for various types of 4d massless staggered fermions along 
(1,1,0). We compare the standard action, a Symanzik improved action called u p6" from 
Ref. H/, and the optimized BA and PD fixed point actions, truncated with mixed periodic 
boundary conditions. 
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